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Small angle neutron scattering (SANS) is an historical technique to study the flux lines lattice 
(FLL) in a superconductor. Structural characteristics of the FLL can be revealed, providing funda- 
mental information for the physics of vortex lattice. However, the spatial resolution is limited and 
all the correlation lengths of order are difficult to extract with precision. We show here that a time 
of flight technique reveals the Bragg peak of the FLL, and also its translational order with a better 
resolution. We discuss the implication of these results for pinning mechanisms in a Niobium sample. 
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INTRODUCTION 

The existence of a flux lines lattice (FLL) in a super- 
conductor of the second kind was initially confirmed by 
Cribier et al using small angle neutron scattering (SANS) 
Afterwards, SANS was an essential technique to 
study the physics of the FLL, including FLL transitions 
@-[5| j pinning mechanisms > moving FLL [H, Q and 
out of equilibrium features [10(. A central question for 
the physics of FLL is the nature of the FLL in the pres- 
ence of unavoidable structural disorder. It is clear that 
the flux lines are forming planes which are ordered at a 
sufficient long range to display Bragg peaks in the diffrac- 
tion pattern. By contrast, a genuine glassy or liquid state 
presents a large spread of Bragg angles and a degener- 
ated orientational order that lead to diffuse rings of scat- 
tering in the reciprocal space. However, it is not clear 
whether the FLL is ordered at very l ong range ll| or is 
fracturing at the intermediate scale [12j. The result is 
important in the field of disordered elastic media where 
novel phases, not existing for the crystalline matter (i.e 
the Bragg Glass) [ijjj], have been proposed. It is also im- 
portant for discriminating between pinning effects, since 
FLL order and bulk critical current can be linked in the 
framework of the elastic theory of collective pinning [3] ■ 
On the other hand, pinning at the surface is also very ef- 
ficient but the associated critical current is not, at least 
not directly, related to the disorder in the bulk FLL struc- 
ture 15j . SANS is the sole technique from which the bulk 



FLL correlation lengths can be in principle extracted |6( . 
Experiments are usually performed with the neutron 
beam applied along the magnetic field, the so-called lon- 
gitudinal geometry. The resolution is good enough in the 
longitudinal direction to extract the longitudinal correla- 
tion length (the straightness of flux lines) with accuracy. 
Note however that this length is directly affected by field 
lines bending due to the demagnetization field [16] or by 



the self field if a transport current is applied [7\ . Its in- 
terpretation is then tricky [Toj]. The scattering intensity 
I(Q) contains the square of the modulus of the structure 
factor S(Q), which is the Fourier transform of the posi- 
tional correlation function C(r). I(Q) is analyzed after 
a radial averaging of the intensity in the detector plane. 
In this direction however, the resolution function of the 
small angle diffractometer is strongly dominating 17j, [18 1 
and a direct information can not be obtained from the 
Bragg peak shape. As a consequence, indications that 
the FLL is in a Bragg glass state arises from the anal- 
ysis of the decrease of the scattered intensity with the 
magnetic field [111 ], and not from the expected power- 
law decay [lj|. For the same reasons, the broadening 
of the peak in the radial direction which gives informa- 
tions on the crystallite size can not be extracted without 
removing the dominating contribution of the resolution 
function. This leads to a large uncertainty in the size of 
crystallites. 

The other geometry is the transverse geometry where the 
beam is perpendicular to the magnetic field, i.e. the high 
resolution direction corresponds to the orientational or- 
der. A major drawback is that the Bragg planes are 
difficult to bring in diffraction conditions for any mis- 
alignment between FLL and the applied field |18| . In 
this geometry, some evidence of a power law decay of 
Bragg peaks are observed in a Niobium of good qual- 
ity, what is consistent with a Bragg Glass phase [19j. 
However, the intrinsic width of these Bragg peaks were 
much larger than the experimental resolution. This im- 
plies that the FLL fractures at an intermediate length 
scale, typically in the micrometer range. These results 
are not expected in the Bragg-Glass scenario of quasi- 
long range order where Bragg peaks have no intrinsic 
width in the regime of power law decay. Complemen- 
tary experimental data should be helpful to clarify the 
situation. 
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We have used another technique to extract the FLL 
crystallite size, using the spectrometer with the usual 
longitudinal geometry but in the time of flight (TOF) 
mode. The time of flight of neutrons between the chop- 
per source and the detector is used to separate in the 
scattering intensity the neutrons of different wavelengths. 
In the conventional steady state mode, a monochromatic 
beam can be selected only with a large wavelength spread 
AA/A ~ 10 — 20 % which has an important contribu- 
tion in the resolution function. In the TOF mode, this 
wavelength spread arising from the pulse and the channel 
widths can be negligible. Another advantage of the TOF 
is that scattering to higher angles 8 which gives a better 
resolution 20] is direct thanks to the large range of wave- 
lengths available in the single experiment. The drawback 
is that the neutron flux is low at large wavelength due its 
the Maxwellian spectrum (the wavelength at maximum 
flux is here A « 2 A). More practical but also interesting, 
it is not necessary to rock the sample through the Bragg 
conditions which are satisfied at a fixed position when A 

is Continuously tuned {Osragg = QBragg^/^ With 8 B ragg 

the scattering angle). The definition of a Bragg peak is fi- 
nally largely improved, i.e. more data points are available 
to describe it what improves significantly the accuracy on 
the peak width. 

SANS data were measured with the PAXY small- 
angle diffractometer (Laboratoire Leon Brillouin, Saclay, 
France) , which can be used in both steady state and TOF 
modes. The sample under study is a slab of pure Niobium 
with critical temperature T c =9.2 K, Ginzburg-Landau 
parameter k « 1 and second critical field i? C 2=0.385 T 
at T=2 K. No peak effect in the critical current was ob- 
served, what is consistent with a sample of the best ho- 
mogeneity (l6| . 

Before the TOF measurements, the FLL was mea- 
sured with the conventional steady state set-up, using a 
mechanical velocity selector which chooses the neutrons 
wavelength A with a Gaussian distribution of full width 
at half maximum (FWHM) AA= 0.11 x A (or a vari- 
ance cr 2 given by AA/2v2/ri2= <j\). The wavelength was 
fixed at A=10 A. I(Q) is obtained after regrouping the 
different rocking curves and after a radial averaging of 
the intensity. We observe in fig.l the Bragg peak asso- 
ciated with the FLL and centered at Qfll — 2tt / (Lfll, 
(d-FLL = 1.155ao, ciq — 1.074-y/(/>o/-E? is the unit cell of 
the hexagonal lattice). The Bragg peak is well fitted by a 
Gaussian because dominated by the spectrometer resolu- 
tion 2l| . This latter is related to the angular distribution 



compared to the one arising from the incident neutron 
beam, a good approximation of the angular variance ag 

is m 



of the incident neutron beam at the detector position, to 
the wavelength distribution AA, and the detector resolu- 
tion AR. Adding the different Gaussian contributions to 
the resolution gives: 



^•2 ~ ~ 2 



Qr 



= (27r/A) 2 ^ + (Q/A) 2 AA 2 /(2V27r^2) 2 . (1) 



The detector resolution contribution being negligible 



sa 1/4(D 1 /2.L 1 ) 2 +1/4[D 2 .(L 1 +L 2 )/2.L 1 .L 2 } 2 

(2) 

with Di, D 2 the aperture diameters, L\ the distance 
between the two apertures and L 2 between D 2 and the 
detector (Here Dl=12 mm, D2=7 mm, Ll=4750 mm, 
L2=6870 mm). 

With our set-up, a Qrea w 6.8 10~ 4 A" 1 at Q = 6.64 
10~ 3 A^ 1 . To measure the intrinsic width of the Bragg 
peak, it is necessary to deconvolve the experimental data 
cjq^ p from <7Q res [I7|, HU. As shown in fig.l, ctq cxp « 
<jQ rcs and we can only conclude that the FLL crystallite 
size has a minimum value of some /^m. We have also per- 
formed full rocking curves of Bragg peaks whose FWHM 
are related to the perfection of flux lines along the field 
direction. We measure resolution limited widths for B= 
2000, 2500, 3000 and 3500 G, indicating a longitudinal 
correlation length of more than 100 fim, i.e. straight flux 
lines. Orientational order was not changing in the field 
range investigated. 

SANS experiments in the TOF mode were then per- 
formed. The velocity selector was removed and a chop- 
per with one slit (12 x 1 mm 2 ) was placed before the 
cryomagnet containing the sample. The chopper to de- 
tector distance was D c h opper =7ll0 mm. The collima- 
tion and the sample to multidetector distance were kept. 
The total pulse width T c /j Opper =450 /is is given by the 
rotation of the 1 mm width slit of the chopper in front 
of the 7 mm diameter sample diaphragm. The flying 
times of neutrons were analyzed in times frame of 256 
TOF channels of r=150 /is each. The scattering inten- 
sity of the sample was recorded at T=2K for different 
magnetic fields in the superconducting state (field cool- 
ing procedure) and a background was measured in the 
normal state (B=4000 G>B c2 ). In the TOF mode, the 
neutrons are recorded as function of the time of flight 
t for different angles 9, then t is converted to effective 
wavelengths, giving 1(6*, A). Each set of raw scattering 
data was corrected for the detector efficiency, the sam- 
ple transmission and the wavelength distribution of the 
incident beam flux by dividing each scattering data by 
I(6>=0, A). Data obtained in the normal state was used 
as the background scattering for data obtained in the 
superconducting state. Typical data showing the FLL 
Bragg peak as function of (6*, A) are shown in fig. 2. In 
TOF mode, the contribution of AA in equation (2) is no 
more coming from the wavelength distribution delivered 
by a velocity selector but arises from the short pulse and 
the TOF channel width. Using A = h/(m n v) (h is the 
Planck constant, m n the neutron mass and v the neutron 
velocity) gives AX chopper = T chopper / (0 .253 x D chopper ) 
and AA T = r/(0.253 x D chopper ) Here, AX chopper and 
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AA r are in A, T chopper and r in /is, D chopper in mm. 
Assuming Gaussian distributions for these two contri- 
butions, the term (AA/A) 2 in equation (1) becomes 
(AX chopper I A) 2 + (AA T /A) 2 . It is much smaller and neg- 
ligible compared to the contribution arising from the an- 
gular divergence of the beam. The position of the peak 
changes with (9, A) according to the Bragg law as shown 
in the inset of fig. 2. The width of the peak decreases 
at large 9, as expected from equation (1). In the TOF 
mode, interestingly, this 9 dependence can be analyzed 
in a single measurement when measuring 1(9, A). 
A typical size of crystallites can be estimated using the 
Scherrer formula 1231: 



A9 



crystal 



K.X/S.cos9 



(3) 



with K the Scherrer constant of the order of unity [24( 
and S the mean size of the crystallite. Note that S is an 
effective length, measured in the direction of the diffrac- 
tion vector. 

Finally, the FLL Bragg peak broadening varies as: 

®Qfll * ^.a e ^ m ) 2 /{d FLL .9) 2 + (2.n.K) 2 / (co S 9.S) 2 

(4) 

Since cos9 ~ 1 at small angle, it can be rewritten in the 
compact form 



A/9 2 + B/S 2 



(5) 



with B a constant of the order of (27r) 2 , and A is function 
of the Bragg planes spacing. 

Despite the gain in resolution offered by the TOF, the 
intrinsic width of the peak still corresponds to a low con- 
tribution compared to the instrumental resolution. The 
equation (5) has in principle two adjustable parameters, 
A and B/S 2 (the instrumental resolution due to the in- 
tercept of the incident beam in the detector plane and 
the crystallite size S). Since the prefactor A in the beam 
resolution contains only geometrical parameters, it can 
be calculated and then fixed at the expected value. In 
such a case, the fit has only one free parameter. In fig. 3, 
we show the comparison between the two fitting proce- 
dures, with one or two adjustable parameters The fits are 
of equivalent quality. From the fitting curves at different 
fields (fig. 4) , we can deduce S as function of the magnetic 
field values. As shown in fig. 5, the order of magnitude 
of S and its field dependence are found similar for the 
two procedures, albeit with larger error bars when using 
the two parameters fit. It certainly reflects the relative 
uncertainty of the extraction of S, but shows also that 
the order of magnitude is correct. In the following, we 
will discuss the case of S deduced from the one param- 
eter fit, but the results and discussion are similar with 
either fitting procedures. 

We observe that S increases with the magnetic field, 
meaning that the FLL becomes more perfect as its den- 
sity increases. S is however slightly larger at 500 G than 



at 1000 G. In a low kappa superconductor such as Nio- 
bium, the first critical field is large, i.e B c \ ~ 1500G at 
2K. The intermediate mixed state where FLL and Meiss- 
ner regions coexist is extended at low field, specially with 
the slab geometry of the sample with a large demagne- 
tizing factor [2|| . The measurements at the lowest field 
are then likely influenced by the increasing presence of 
Meissner domains in the sample. 

In the mixed state, the critical current density J c (bulk) 
is related to the positional correlation length R c in the 
Larkin Ovchinikov (LO) model which is valid for 
short lenght scales displacements, typically r p , the one 
of the pinning potential [26|. It has been pointed out 
that a diffraction experiment probes scale of the order 
of the lattice spacing [27[. The correlation length S 
can not be directly compared with R c , but should be 
rescaled as S 



Rc{d F LL/r p ) 2 



27|. 



is expected to 



vary between the flux line core size £ for isolated flux 
lines up to some fraction of the lattice spacing <Ifll 
when flux lines largely overlap (jjj . We have calculated 
J c (bulk) sa C e6 r p /(4:BRl) [l4[ with the two limiting val- 
ues of r p . The shear modulus Cqq formula was given in 
29], and £(2K) w ((f> /2irB c2 ) 1 / 2 « 29 nm. J c (bulk) is 
finally reported in Fig. 6, for r p = £ and for r p = d F LL- 
The next step is to compare with the real critical current 
J c (exp) of the sample. It was measured on a small piece 
of the same Nb sample (width x length x thickness = 
0.1 x 0.3 x 0.02 cm 3 ), using the irreversible part of the 
magnetization and applying the Bean model for a slab ge- 
ometry. The resulting J c (exp) is shown in Fig. 6. Clearly, 
J c (bulk) is very different than J c (exp) for B/B c2 > 0.4, 
in the regime where our Niobium sample is clearly in the 
pure mixed state, whatever the r p value. 

We conclude that a large amount of critical current is 
not coming from a bulk pinning contribution of the LO 
type. Another possibilty is that we are measuring a crys- 
tallite size unrelated to any bulk pinning mechanism, for 
reasons which are not clear. In soft superconductors of 
the second kind, the bulk contribution to the pinning can 
be very low and most of the critical current can arise from 
a surface origin. In particular, surface pinning is known 
to be a realistic source of pinning in Niobium [lil l30j . 
This pinning mechanism is based on equilibrium equa- 
tions and boundary conditions for flux lines over a realis- 
tic (rough) surface pSlj . The surface contribution to the 
pinning gives a critical current i c (A/m) = e.sind where e 
is the vortex potential (i.e. the reversible magnetization) 
and 9 is a critical angle characterizing the surface rough- 
ness. i c is a superficial current but can be rewritten as 
a critical current density J c (surf)(A/m 2 ) = 2i c /t where 
t is the sample thickness and the factor 2 stands for the 
two surfaces perpendicular to the magnetic field, stan- 
dard values of 9 are around few degrees. To estimate i c , 
e was computed as function of B using numerical calcu- 
lation following 32j. For this calculation, we have used 
a Ginzburg Landau parameter k = 1 and B c2 = 0.385 
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T as observed experimentally. The result is reported in 
Fig. 7. A direct observation is that the magnetic field 
dependence of e and J c (exp) are quite similar, as ex- 
pected if the critical angle is not much depending on the 
value of the field. Fitting the experimental J c (exp) of 
our Niobium sample with the surface pinning expression 
J c (surf) = 2/t x e x sind gives a very good agreement 
with a critical angle 9 — 4 deg (Fig. 8). This is a large, 
but still reasonable value for a bulk sample with unpol- 
ished surfaces [IH . We do not exclude that edge currents 
play also a role for the critical current. Finally, the good 
agreement between the surface pinning model and the 
experimental data, in addition to the neutron scattering 
experiment, allow to conclude that surface critical cur- 
rents are likely dominant in this sample. 
Another important result of the TOF experiment is that 
the FLL is fracturing in the micrometer scale, as found 
using reverse monte carlo analysis in a different geometry 
[l^ | . This result is a priori not expected in the framework 
of purely elastic models where quasi long range ordering 
is expected in clean samples such as Niobium. Note that 
we have used field cooled, zero field cooled and the so- 
called shaking procedures to induce a better positional 
order of the FLL [34j]. The widths of the peak are not 
changed (within resolution), meaning that unpaired dis- 
locations are likely not responsible for the Bragg peak 
width and that the FLL is close to equilibrium. 

To conclude, we have measured the crystallite size of 
the FLL in Niobium using neutron scattering with a time 
of flight mode. Crystallite sizes are in the micrometer 
range, and increase with the field. These results show 
that the FLL positional order gets better when the flux 
line density increases. The crystallite sizes are found 
unrelated to the critical current using a bulk collective 
pinning approach. This implies other source of pinning 
unrelated to the bulk FLL structure, likely surface pin- 
ning which is found to describe quantitatively the critical 
current values. 
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FIG. 1: Normalized intensity as function of the scattering 
vector Q (A = 10A, T=2K)). The first Bragg peak is fitted by 
a Gaussian of center Q= 0.0066 A" 1 and FWHM 0.00165 ± 
0.00004 A -1 (solid line), very close to the calculated FWHM 
due to the resolution (0.00160 A" 1 ). 




FIG. 2: Intensity of the FLL Bragg peak / as function of 
^neutrons (time of flight mode) . Each Bragg peak corresponds 
to a different (0=0.201, 0.229, 0.249, 0.271, 0.291, 0.312, 
0.333, 0.354, 0.375, 0.396, 0.412, 0.437 deg). In the inset is 
shown Xp e ak, the center of the Bragg peak, as function of the 
angle 0. The solid line is the Bragg law, with QBra gg ~ 7.47 
10~ 3 ± 0.02 A" 1 (this is the value expected for the hexagonal 
lattice at 2500 G within the resolution of the magnet). 
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FIG. 3: uq f ll as function of the angle of scattering 9 (T=2K, 
B=2000 G). In both graphs, the solid line is a fit with the 
equation <7q ~ A/6 2 + B/S 2 , with one free parameter 
B/S 2 (left) or with two free parameters A and B/S 2 (right) 
(see text). R 2 is the coefficient of determination. 




9 (deg) 



FIG. 4: &q f ll as function of the angle of scattering 8 (T=2K, 
different field values) . The solid line is a fit with the equation 
oq fll ~ A/9 2 + B/S 2 , with B/S 2 as a free parameter(see 
text). For clarity, the fit with two parameters is not shown. 
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FIG. 5: Size of the vortex crystallite as function of the reduced 
magnetic field B/B c2 (T=2K, B c2 =3850 G) obtained with the 
single parameter fit or with the two parameters fit. 
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FIG. 6: Experimental critical current measured at T=2K 
(stars), compared with the critical current calculated with the 
LO model and S values as explained in the text (empty points 
correspond to r p = £ and empty triangles to r p — (Lfll)- 
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FIG. 7: Potential e (or reversible magnetization) computed 
using the Brandt approach [32| as function of reduced mag- 
netic field in a semi -log scale (k—1 and B C 2 = 0.385 T). In the 
inset is shown the same graph in a linear scale. The dotted 
line is the Abrikosov line. 
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FIG. 8: Experimental critical current measured at T=2K 
(stars), compared with the critical current calculated with 
the surface pinning model with a critical angle 8 — 4 deg 
(plain points). 



